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Abstract 

During heterogeneous nucleation within a metastable mushy-zone, several geometries for 
nucleation site development must be considered. Traditional spherical dual cap and crevice 
models are compared to a truncated dual cap to determine the activation energy and critical 
cluster growth kinetics in ternary Fe-Cr-Ni steel alloys. Results of activation energy results 
indicate that nucleation is more probable at grain boundaries within the solid than at the solid- 
liquid interface. 


Introduction 


Undercooled hypoeutectic ternary Fe-Cr-Ni steel alloys solidify in a two-step process known as 
double recalescence [1]. The metastable ferritic bee-phase forms first with subsequent 
transformation to the stable austenitic fee-phase. The delay between these two nucleation events 
represents the incubation time for the stable phase and this delay is a strong function of melt 
convection. In classical nucleation theory [12-14], the steady-state nucleation rate, Is, and time 
dependent nucleation rate, I, are related by the equations 
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where x is a characteristic incubation time, t the observed delay time, Io a pre-exponential factor, 
AG* the Gibbs free energy for formation of a critical nucleus of n* atoms of a pure material, and 
T the transformation temperature. The Boltzmann constant has a value of k B = 1.38 x 1 0 
J/atomK. Steady-state nucleation does not become appreciable until AG* > 60 k B T while 
transient nucleation does not become appreciable until t » x. Turnbull estimated that for 
condensed phases the incubation time was a function of the size of the cluster and the rate at 
which atoms cross the interface, P*, between the surrounding matrix and the cluster [2], 
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The attachment rate [3* has variously been related to the jump frequency [2] and the lattice 
diffusivity [3-5]. By invoking the principle of time reversal, where the statistical fluctuations in 
cluster size follow the same path during growth and decomposition, Feder et al [6] determined 
the incubation time to be 
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Russell [7] evaluated condensed phase nucleation for binary systems to account for the influence 
of solute partitioning and found that clusters approaching the critical size were surrounded by an 
enriched solute shell, contrary to expectation, and that diffusivity of the slower moving species 
controls the attachment rate with an interchange frequency P'. 
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where x is the number of atoms jumping a jump distance a 0 and the subscripts on the diffusivity 
D represent solute A and solvent B. The number of atoms jumping can be evaluated from the 
surface concentration C A - the interface area As, and the area per atom. A key finding for multi- 
component systems is that the rate controlling step for linked flux evaluations is replacement of 
shell atoms of the controlling species and not the interfacial jump frequency. 

In an effort to identify the mechanism for stable-phase nucleation in ternary Fe-Cr-Ni steel 
alloys, Koseki [1] showed that a dual cap cluster forming within the metastable solid along a 
bcc-bcc boundary was more probable than (1) homogeneous nucleation of a sphere in the melt, 
(2) heterogeneous nucleation as a single spherical cap protruding from the bcc-solid and into the 
melt, (3) heterogeneous nucleation of a spherical cap protruding into the solid, and (4) 
heterogeneous nucleation of a dual cap along the solid-liquid interface. In these analyses, the 
time to form a potential nucleation site is taken as zero. Homogeneous nucleation in the melt 
requires no site development time. Homogeneous nucleation within the solid and heterogeneous 
nucleation in the form of a spherical cap at the solid-liquid interface requires only that the pre- 
existing solid exist. Since the solid/liquid interface forms immediately during primary 
recalescence there is no delay for site formation. 

A transient site, such as the solid-solid grain boundary dual cap, is one that forms sometime after 
primary solidification. Two processes become important in the evaluation of transient site 
development. First, we must assess how long it takes to form a suitable heterogeneous nucleation 
site and then, second, we must define how long it takes for a critical nucleus to grow after a 
suitable site has formed. The sum of these two events determines the incubation period which 
corresponds to the delay time between primary metastable recalescence and subsequent initiation 
of the conversion to the stable phase. 

Hanlon [8] proposed that fluid drag and dendrite bending can lead to secondary arm collision 
which will create the required bcc-bcc grain boundary. Matson [9] extended these analysis by 
showing that nucleation of the stable phase could only occur once a suitable site had formed on, 
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or in, the pre-existing metastable solid thus explaining how convection can influence predictions 
using classical nucleation theory. Sites that form due to collision include the bcc-bcc grain 
boundary and the liquid-filled crevice formed by the intersection of two parallel cylindrical arms. 
Since the predicted delay time required for cluster growth for both the dual cap and crevice 
clusters was several orders of magnitude faster than the observed delay, the formation of the 
grain boundary crevice, and not cluster growth, was identified as the controlling process. 

One additional geometry needs to be considered. The current paper looks at the potential for 
interaction between a dual cap and the liquid due to growth near the grain boundary/liquid 
interface. Dual cap cluster formation was primarily limited due to diffusion along the grain 
boundary since diffusion through the solid was sluggish and thus grain boundary diffusion 
dominated. The activation energy was low due to the favorable low surface/volume ratio. 
Crevice cluster formation was significantly more rapid due to enhanced diffusion through the 
surrounding liquid but the activation energy was higher due to an increase in surface/volume. 
The intersection of a dual cap with the liquid surface combines these processes and must be 
evaluated to complete the picture. The geometry that forms is known as a truncated dual cap. 

Geometry of a truncated dual cap 

A dual cap consists of two identical spherical sections with a characteristic intersection angle, 29, 
defined by the relative interfacial free energies, y. Along the grain boundary mirror plane, 
equilibrium is maintained if 

7 mm = 2 7 ms cos 0 (5) 

with subscripts MM for metastable/metastable and MS for metastable/stable. The activation 
barrier at the critical radius, AG*, is evaluated based on equating surface and volume terms, A 
and V respectively, at the critical cluster radius r* where dG /d n goes to zero such that 

AG *L. = 2 A MS 7 ms - A MM 7 mm ~ V AG V (6) 

AH AT 


with AH the enthalpy of fusion, AT the undercooling set by the melting points of the metastable 
and stable phases, and T m the melting point of the stable phase. The area and volume terms are 
readily evaluated from equations for spherical sections. If we assume typical values for 
thermophysical properties for steel alloys[9] to be y M M = 0.468 J/m 2 , y M L = 0.212 J/m 2 , 
y SL = 0.302 J/m 2 , AG V = 2.56 xlO 7 J/mol, D GB = 1 x 10" 10 m/s 2 , D s = 1 x 10" 14 m/s 2 , 
Dl = 3 x 10" m/s", and a grain boundary thickness of t = 5 A, then the angle is 29 = 41.2°, the 
activation barrier is ag *4t =102 for a nucleus containing 15000 atoms with a delay time of 
x = 5.9 xlO" 8 sec. 

Note that the value for the interfacial free energy between stable and metastable phases is highly 
speculative given the inability to experimentally define this value. Thus a value of 
Yms = 0.250 J/m 2 was selected as representative. Throughout the remainder of this paper, the 
errors associated with this selection are investigated. 
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When a dual cap intersects the liquid, the common circular basal surface is truncated. This is 
shown in Figure 1 where the overall shape of the transformed volume is seen in (a) and each of 
the three sub-volumes which comprise the shape are visually blown apart in (b). Note that top 
and bottom represent the converted section of the dual cap within the metastable solid while the 
pointed ellipsoid extends out into the liquid. 



Figure 1. Geometry of a truncated dual cap. (a) overall volume (b) exploded view to show 
upper and lower truncated cap sections and the ellipsoid protruding out into the 
liquid (c) truncated cap for typical properties and yms = 0.25 J/m 2 (d) truncated 
cap for typical properties and yms = 0.37 J/m 2 

In Figure 1 the influence of the value of the interfacial free energy between metastable and stable 
phases is illustrated where in (c) the interfacial free energy is set at Yms = 0.25 J/m" while in 
(d) Yms = 0.37 J/m 2 . Of particular interest is the strong influence of the value of the y M s 
parameter on the ellipsoid dimensions. As seen in Figure 1(c), at low values the interfacial 
properties dominate and the point where metastable solid/metastable solid/liquid are in 
equilibrium is pronounced. In Figure 1(d) the volume properties dominate and the point is less 
pronounced. This point is of critical concern because at this location surface frustration, a local 
instability caused by anisotropy of surface tension, causes modeling of the geometry to break 
down. Isotropic interfacial energies require that the boundary between metastable solid and 
liquid have the same wetting angle which means that the “ellipsoid” should, in fact, be a 
spherical cap section with a-axis and b-axis the same. This is certainly not true in (c) but the 
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geometry comes closer to this shape in (d). Surface frustration is due not only to the elliptical 
cross-section of this face but also due to the definition of surface angles. 

Along the solid/solid interface (left side of the shape) equation 5 is used to define the interior 
angle 29 relative to the horizontal mirror plane as before. At the point where the metastable 
solid-solid boundary intersects the liquid-stable solid boundary (the point in the ellipsoid) the 
interfacial free energy balance results in the definition of two equations in two unknowns along 
horizontal and vertical axes 

YlM = r MS COS 9 S + 7 IS cos e L ( 7 ) 

Ym S sin d S = TlS sin °L (8) 

These define the angle the truncated section makes to the vertical, 9s, and the angle the exterior 
surface of the ellipsoid makes to the vertical, 9l ■ Note that ideally the angle should be identical 
along the a-axis and b-axis of the ellipsoid but surface frustration makes this not possible due to 
the required existence of the pointed intersection. To estimate the surface area and volume it 
was assumed that the shape could be represented by rotating the destroyed metastable solid- 
liquid interface about the solid-solid line resulting in a shape similar to half a football. 

The truncated dual cap and ellipsoid shapes cannot be characterized readily analytically and a 
graphical approach was used to estimate areas and volumes. Using the predicted characteristic 
lengths and angles, a model was developed in SolidWorks® to evaluate areas and volumes as a 
function of the interfacial free energy Yms- Defining the base of the truncated cap as A MM , the 
front as A M l, the surface as A M s, and the exterior of the ellipsoid as A S l, the activation barrier 
can be evaluated from 

AG *L = 2A msYms + A slYsl - 2 A mlYml - A mm Ymm ~ V AG v (9) 

Using typical thermophysical properties presented earlier, 9 = 41.2°, 9s = 81.2°, 9l = 54.9°, and 
the ellipsoid has major axis length a = 19.9 x 10 9 m and minor axis length b = 3.3 xlO 9 m. The 
activation barrier is calculated to be AG */ kT = 3315 for a nucleus containing 156000 atoms with a 
delay time of x = 2.8 xlO 9 sec. The results of calculations of this type for all three mechanisms 
as defined over the entire range of possible interfacial free energies are presented in Figure 2 and 
Figure 3. 


Discussion 

Prediction results show that all three mechanisms are possible but only for low values of the 
interfacial free energy. At these energies, the shape of both the dual cap and truncated dual cap 
are thin saucer-shaped volumes due to the small interior angle along the metastable solid grain 
boundary. The truncated portion thus has a highly anisotropic elliptical shape with major axis to 
minor axis ratio varying from 11 to around 6 over the range of interfacial free energies 
investigated. The sensitivity to selection of values for interfacial energies can readily be 
investigated beyond what has been presented here as part of a larger evaluation of error 
associated with use of this type of technique. For example, if the value of the interfacial free 
energy of the metastable solid-solid boundary is reduced from the nominal value of y M M = 0.468 


197 



2 2 

J/m to Ymm = 0.212 J/m" (the value for y M l) then major to minor axis ratio is reduced to around 
1.6 and is insensitive to the selection of Yms- This ratio represents geometries closer to the 
desired isotropic, circular shape. Note that all interfacial free energy values have a high inherent 
uncertainty and thus the results should be viewed as indicating trends and not as numerical 
predictions. Results are ordinal and not cardinal. 
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Figure 2. Energy barrier for nucleation of each cluster shape over the range of possible 
interfacial free energy values. 


Figure 2 also clearly shows that the truncated cap geometry does successfully combine elements 
of both the solid-solid dual cap mechanism and the solid-liquid crevice mechanism. Predictions 
for the truncated dual cap tend to lie between the extremes calculated for the other two 
mechanisms however if all were operational this mechanism would dominate since the fastest 
possible mechanism is favored and the truncated cap is kinetically fastest over the entire range, 
as seen in Figure 3. However, since the energy barrier is higher for the truncated dual cap than 
for the embedded dual cap it is unlikely to be operative except at interfacial free energies less 
than a value of around Yms = 0-24 J/m . 

One additional weakness of the approach used to estimate the surface area and volume for the 
pointed ellipsoid section for this geometry is that in order to obtain a tractable solution for the 
ellipsoid form, the foot-ball shaped rotation element was selected. This technique was used 
because SolidWorks® can readily generate this shape. Unfortunately, this ignores the calculated 
angle 0 L = 54.9° for Yms = 0.25 J/m which should be constant along the three-phase interface. 
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Figure 3. Kinetics of cluster formation for each shape over the range of possible interfacial 
free energy values. 

Note that the rotated shape has a continuously changing surface angle. The extremes occur at the 
point along the major axis where the angle is 9 [ nun = 18.9°, and at the top along the minor axis 
where the angle is 0 [ max = 90°. The effect of this change in angle is a continuously changing 
curvature of the surface which, as previously noted, results in surface frustration. This warping 
of the ideal geometry also results in a slight overestimation of the surface area and volume - 
estimated to be on the order of 37% and 58% respectively for these conditions - with the error 
decreasing because the difference between major and minor axis decreases significantly with 
increasing y M s- One interesting result of this surface frustration is that the decreased angle at the 
point strains the equilibrium at the three-phase junction and creates an additional resistance to 
cluster formation. This effect is commonly observed by children blowing up a balloon - if a 
strain is induced along one axis by pinching, the increase in anisotropic surface tension only 
allows expansion of the balloon along orthogonal axes. Thus for our application there is a strain 
term which must be added to the free energy balance and growth along the major axis is 
hindered. The polarity of this term suggests that the resulting activation energy will be greater 
than anticipated. These observations further strengthen the conclusion that clusters will form 
along the grain boundary within the metastable solid and not at the intersection of the metastable 
boundary and the undercooled liquid. 

These theories must be validated over a wide range of convection conditions as is only possible 
through electromagnetic levitation experimentation in microgravity [10]. 
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Summary 


Evaluation of the predicted delay times for all mechanisms shows that the incubation time to 
grow a critical cluster is several magnitudes faster than observed. This supports the theory that, 
following formation of the metastable array, the controlling kinetic mechanism is formation of a 
site where nucleation may occur rather than growth of a critical cluster. Evaluation of the energy 
barrier indicates that over a select and limited range of interfacial free energy values, all of the 
mechanisms - spherical cap, crevice nucleation, and truncated cap - are possible. However, 
since the barrier for formation of a truncated cap is greater than that for a dual cap for all 
interfacial energies evaluated in this study, cluster formation is more likely along a grain 
boundary within the solid rather than at the intersection of the boundary and the liquid. 
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